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ABSTRACT 



IC 2118, also known as the Witch Head Nebula, is a wispy, roughly cometary, ~5 

degree long reflection nebula, and is thought to be a site of triggered star formation. In 

l/y order to search for new young stellar objects (YSOs), we have observed this region in 

7 mid- and far-infrared bands using the Spitzer Space Telescope and in 4 bands in the 

M-r optical using the U. S. Naval Observatory 40-inch telescope. We find infrared excesses 

in 4 of the 6 previously-known T Tauri stars in our combined infrared maps, and we 
V— i ' 
-*— > ■ find 6 entirely new candidate YSOs, one of which may be an edge-on disk. Most of the 

YSOs seen in the infrared are Class II objects, and they are all in the "head" of the 

nebula, within the most massive molecular cloud of the region. 



Subject headings: stars: formation - stars: circumstellar matter - stars: pre-main 
sequence - ISM: clouds - ISM: individual (IC 2118) - infrared: stars - infrared: ISM 



^1- 

q 

Spitzer Science Center/Caltech, M/S 220-6, 1200 E. California Blvd., Pasadena, CA 91125 

2 Current addres: European Southern Observatory, Alonso de Cordova 3107, Casilla 19001, Vitacura, Santiago 19, 
Chile e-mail: sguieu@eso.org 

3 U.S. Naval Observatory, Flagstaff Station, 10391 West Naval Observatory Rd., Flagstaff, AZ 86001-8521 

4 Oil City Area Senior High School, Oil City, PA 



5 Bassick High School, Bridgeport, CT; currently New Jersey Astronomy Center, Raritan Valley Community 
College, Somerville, NJ 

6 Lincoln High School, Stockton, CA 

7 Luther Burbank High School, San Antonio, TX 

8 McDonogh School, Owings Mills, MD 

'Pomona College, CA 

10 NASA Exoplanet Science Institute, IPAC, M/S 1002-22, 1200 E. California Blvd., Pasadena, CA 91125 



Introduction 



Our general understanding of how and why it is stars form rests upon the answers to a set 
of more basic questions: given a particular locale, how, why, and when did stars form there? 
Current theory tells us that in some regions, stars simply collapse under the influence of their own 
natal clouds' gravity, but in others, the process may be heavily influenced by external events such 
as nearby supernovae, which may prompt parts of the cloud to form stars. Many people have 
searched for evidence of triggered star formation on large and small scales (e.g., Briceho et al. 
la and references therein). In the cases where triggered star formation is suspected, a complete 
inventory of the young stars formed by the cloud is useful specifically to (a) search for evidence of 
sequential star formation within the cloud and with respect to other nearby clouds, and (b) within 
the individual cloud, understand the initial mass function (IMF), the total mass in the cloud vs. 
stars, and the star formation efficiency The former strengthens (or weakens) the case for having 
found legitimate triggered star formation, and the latter allows closer study of the initial conditions 
of star formation (e.g., Ballesteros-Paredes et al.[l|). If the specific cloud does turn out to be forming 
stars as a result of triggering, then a detailed comparison of the IMF and star formation efficiency 
for a triggered cloud versus a cloud that has collapsed under its own gravitational influence will 
shed further light on the mechanics of star formation itself. 

Finding all of the young stars in any given cloud in the optical can be problematic due to 
extinction from the cloud itself. However, the Spitzer Space Telescope ()50|) provides an excellent 
platform for surveying star- forming regions in the mid-IR and far-IR, enabling stars with IR excesses 
to be identified in a relatively straightforward fashion, even those embedded fairly deeply in the 
cloud. Because Spitzer efficiently maps large regions of sky, it is a particularly useful tool for 
surveying clouds subtending relatively large areas. 



IC 2118 (J22| ;|24|). the Witch Head Nebula, is ~5 degrees long; its "wind-blown" appearance is 
similar in the optical and infrared - see Figure [1] for the wide- field optical image from the Palomar 
Observatory Sky Survey. It has been cited as an example of triggered star formation. The specific 
distance to the obj ect is somewhat controversial. If it is close to ~400 pc, the trigger could be 
the Trapezium (126l ). If, instead, it is closer, some ~200 pc away, the possible trigger is the Orion- 
Eridanus superbubble (1221: \2m ) . and the cloud is likely to be additionally excited by /3 Orionis (Rigel, 
at 264 pc; van Leeuwen 143) . It is possible that interactions with both the expanding superbubble 
and the winds from Rigel (both to the East of IC 2118) create the cometary shape of the nebula. 
For the rest of this paper, we take the distance to IC 2118 to be 210 pc (|22j), although we comment 
more on this distance uncertainty below. 



Some previous authors have noted star formation in this region; (122 ) summarizes all previously 
identified objects here, including red nebulous objects found by (J9|). With a 12 CO survey of the 
IC 2118 region, ([22) detected six molecular clouds which include MBM 21 and MBM 22 (G 208.4- 
28.3 and G 208.1-27.5, respectively). Kun and collaborators ()22j; [23j; |24|) conducted an objective 
prism survey with spectroscopic follow-up to search for new T Tauri stars in the IC 2118 region. 



They concluded that there are 11 young stars in the vicinity, 8 of which they discovered, and 6 of 
which are located within (projected onto) the higher-density nebulosity of the cloud; see Table [Q 



(where all the types come from Kun et al. (|24l )) and Figure [TJ Based on this census of YSOs in 



IC 2118, the star formation efficiency here was taken to be ~3%. 

We used Spitzer to survey the heart of the IC 2118 region, with the Infrared Array Camera 
(IRAC; Fazio et al. 1 14 ) at 3.6, 4.5, 5.8, and 8 microns, and the Multiband Imaging Photometer for 
Spitzer (MIPS; Rieke et al. |4Cj) at 24, 70, and 160 microns. In this paper, we use these data to 



search for candidate young stellar objects (YSOs) with infrared excesses. We use our additional 



optical and Two- Micron All-Sky Survey (2MASS; Skrutskie et al. I43l) photometric data in an effort 
to segregate YSOs from background galaxies. While one previously-known T Tauri star appears 
without an excess in our MIPS map alone, we find infrared excesses in 4 of the 5 previously-known 
T Tauri stars in our combined IRAC & MIPS maps, and we find 6 entirely new candidate YSOs. 

The observations and data reduction are described in Sj2j We select YSO candidates using 
Spitzer colors in £[3l and discuss their overall properties in ^J Finally, we discuss some wider 
implications and summarize our main points in [j5j 



2. Observations and data reduction 

2.1. Spitzer Observations 

We observed IC 2118 with Spitzer using both IRAC and MIPS. These observations were 
obtained as director's discretionary time (DDT) as part of the Spitzer Observing Program for 
Teachers and Students. Most of the data were part of programs 235 and 266 (PI: T. Spuck), 
obtained in Spring 2005 and 2006, respectively; additional observations of small portions of the 
map were obtained as part of program 462 (PI: L. Rebull) in Fall 2007 and Spring 2008 in an 
effort to increase the legacy value of the data set. The IRAC observation was broken into nine 
astronomical observation requests (AORs), and the MIPS observation into three; the associated 
AORKEYq^l are given in Table [2j The IRAC observations used the 12-second high-dynamic-range 
(HDR) mode such that each exposure consists of a 0.6 sec and a 12 sec frame. There were three 
dithers per position to minimize instrumental effects, stepping by 260" for each portion of the map. 
The IRAC observation covers ~1.6 deg 2 in all four IRAC bands (3.6, 4.5, 5.8, and 8 /im); Figure [2] 
shows a three-color IRAC mosaic made using 4.5, 5.8 and 8 //m. The outlines of our MIPS and 
optical observations are included for reference. 

The MIPS observations were fast scans, and each scan leg is offset by 160" (half an array). 
There is complete coverage at 70 /im and two scan legs per position at 24 /im. We did not anticipate 



1 AORKEYs are the unique 8-digit identifier for the AOR, which can be used to retrieve these data from the 
Spitzer Archive. 
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Fig. 1. — Palomar Observatory Sky Survey (POSS) red image of IC 2118, with the eleven previously- 
identified T Tauri stars indicated (red circles), as well as the six new candidate YSOs found here 
(blue x). The outlines of our IRAC (red), MIPS (black), and optical (green) observations are 
included for reference. 



having viable 160 /im data, so we did not plan on complete coverage at this bandpass. However, 
the observations were conducted in cold MIPS campaigns, which results in viable data even if the 
coverage is incomplete. The MIPS observation covers ~4 deg 2 ; the 24, 70, and 160 |im full mosaics 
are shown in Figure [3) 



2.2. IRAC Data Reduction 

For the IRAC data, we started with the Spitzer Science Center (SSC) pipeline-processed basic 
calibrated data (BCDs), version S18.7 for program 235 and S18.5 for the two other programs; for 
our purposes, there is no difference between these pipeline versions. We further processed the BCDs 
with the IRAC artifact mitigation software written by S. Carey and available on the SSC web site. 
Mosaics were constructed from the corrected BCDs using the SSC mosaicking and point source 



extraction software package MOPEX ((32j). The mosaics were computed to have a pixel scale of 
1".22 pixer 1 . 

We ran the APEX portion of MOPEX for source detection on the mosaics; we performed 
aperture photometry on the combined long- and short-exposure mosaics separately using the APEX 
output and the "aper" IDL procedure from the IDLASTRO library. We used a 2 pixel aperture 
radius and a sky annulus of 2-6 pixels. The (multiplicative) aperture corrections we used follow the 
values given in the IRAC Data Handbook: 1.213, 1.234, 1.379, and 1.584 for IRAC channels 1, 2, 
3, and 4, respectively. Flux densities have been converted to magnitudes with the zero magnitude 
flux densities of 280.9 ± 4.1, 179.7 ± 2.6, 115.0 ± 1.7, and 64.1 ± 0.9 Jy for channels 1, 2, 3, and 
4 respectively, as given in the IRAC Data Handbook. 

To make a realistic estimate of the internal uncertainties, we obtained photometry from the in- 
dividual BCD images for the candidate members. We found that the scatter within that distribution 
is generally of order 1%. Absolute uncertainties, arising from source variability and uncertainties 
in the instrument calibration (primary calibrators, aperture corrections, etc.), will be an additional 
error. This absolute calibration uncertainty is estimated to be at least ~1.4%. For the longer 
integration times, we estimated average photometric errors in the four IRAC channels to be 0.02 
mag for bright stars, increasing to 0.1 mag at 17.2, 16.6, 14.4, and 13.8 mag for channels 1 to 4 
respectively. 

The APEX source detection algorithm has a tendency to identify multiple sources within the 
PSF of a single bright source, which can cause significant confusion at the bandmerging stage. 
Because our flux densities are determined from aperture photometry with a 2 pixel aperture radius 
that is not deblended, any object with a companion within 2 pixels is not actually resolved. The 
photometry lists were therefore purged of multiple sources prior to bandmerging; objects with an 
accidental companion within 2 pixels had the lower signal-to-noise source removed. 

We extracted photometry from the long and short exposure mosaics separately for each channel, 
and merged these source lists together by position using a search radius of 2 pixels (2". 44) to obtain 
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Fig. 2. — Three-color view of IC 2118 with 3.6 //m in blue, 4.5 /im in green and 8 fim in red. 
Outlines of our optical (black line) and MIPS (green line) surveys are also plotted. 
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Fig. 3. — The 24 /mi, 70 //m and 160 /mi MIPS mosaics of IC 2118. Outlines of our optical (black 
line) and IRAC (red line) surveys are also included on the 24 /mi mosaic. The arcs are from circles 
that are centered on the Trapezium (green circle, 7.35° radius) and on Rigel (blue circle, 4.48° 
radius). 



a catalog for each channel. The magnitude cutoff transition, where we used the long exposure 
rather than the short exposure photometry, corresponds to magnitudes of 11, 10, 8.4, and 7.5 
mag for channels 1 through 4, respectively. We obtained ~40,000 sources in channel 1, ~34,000 
sources in channel 2, ~7500 in channel 3, and ~5000 sources in channel 4. These relative numbers 
are consistent with expectations from other star-forming regions and with the sensitivities of the 
channels. 



2.3. MIPS Data Reduction 

For MIPS, we started with the SSC-pipeline-produced BCDs, version S16.1. The 24 /mi data 
were processed differently from the 70 and 160 /jm data. 

For 24 /im, we used MOPEX to construct a mosaic from the BCDs, with the pixel scale set at 
2.45", very close to native pixel size. We used the APEX 1-Frame routines and the point-response 
function (PRF) provided on the SSC website to perform point-source PRF-fitting photometry on 
the mosaics. We determined via visual inspection that 24 /mi sources with an APEX-extracted 
signal-to-noise ratio (SNR) below 7 were not reliable; therefore, we did not include those sources in 
the following analysis. Moreover, we inspected individual sources and removed 13 more clearly false 
detections (arising from nebulosity or diffraction patterns around bright stars). This resulted in a 
catalog of 1082 point sources ranging from 0.9 mJy to 700 mJy. The zero point used to convert flux 
densities to magnitudes is 7.14 Jy, based on the extrapolation from the Vega spectrum as published 
in the MIPS Data Handbook. We took the error from the signal-to-noise ratio returned by APEX 
as the baseline statistical error. However, the true uncertainty is affected by flat fielding issues, 
etc., and is estimated to be ~4%. We added a 4% flux error in quadrature to the APEX-derived 
error to estimate our final 24 /mi uncertainties. 

At 70 and 160 /zm, the individual unfiltered frames (more than 5400 BCD files) were combined 
using MOPEX, setting their pixel scales to their native values of 9.8" and 16" for 70 and 160 /mi, 
respectively. The final mosaic at 70 /im has a peak surface brightness of ~ 100 MJy sr" 1 , which 



is within the linear regime of the 70 /xm array (Il6l ): for the 160 /mi mosaic, the value is ~ 200 



MJy sr" 1 which is also in the linear regime of the array (44J), and therefore we expect at most a 
20% absolute flux density uncertainty. The data were further processed to maintain the structure 
of the extended emission and preserve its calibration. For the 160 /mi mosaic, we applied a two- 
dimensional (3x3 native pixel) median boxcar filter to interpolate across missing ("NaNs" or "not 
a number") data due to the incomplete coverage (see above). This step affects only those pixels 
near the "NaNs" by a few percent. For both Ge channels, we destriped the image using a ridgelet 
algorithm (J. Ingalls, private communication), which conserves flux and reduces the noise. 

For the 70 /mi point sources, we performed point-source photometry in the same fashion as we 
did for 24 /mi, using APEX 1-Frame on the mosaic constructed as described, and the SSC-provided 
PRF. Sources with SNR below 8 were rejected, and we carefully removed obvious false detections 



Table 1. Previously-identified T Tauri stars in IC 2118 



RA (J2000) Dec (J2000) 



spectral typc a 



other names 



IRAS 04591-0856 


05 01 30.2 


-08 52 14 




2MASS 05020630-0850467 


05 02 06.3 


-08 50 47 


M2 IV 


RXJ 0502.4-0744 b 


05 02 20.8 


-07 44 10 




2MASS 05060574-0646151° 


05 06 05.7 


-06 46 15 


G8: 


2MASS 05065349-0617123 c ' d 


05 06 53.5 


-06 17 12 


K7 IV 


2MASS 05071157-0615098 c ' d 


05 07 11.6 


-06 15 10 


M2 IV 


2MASS 05073016-0610158°' d 


05 07 30.2 


-06 10 16 


K6 IV C 


2MASS 05073060-0610597 c ' d 


05 07 30.6 


-06 11 00 


K7 IV C 


RXJ 0507.8-0931 


05 07 48.3 


-09 31 43 




2MASS 05094864-0906065 


05 09 48.6 


-09 06 07 


G8 


2MASS 05112460-0818320 


05 11 24.6 


-08 18 32 


MO 



HHL 17, G13 

2MASS 05022084-0744099 

(may not be a member of IC 2118; see Kun ct al. 24) 

IRAS F05047-0618 

IRAS 05050-0614, Kiso Ha A0974-19, RNO 37 

Kiso Ha A0974-20, RNO 37 

2MASS 05074833-0931432 

(may not be a member of IC 2118; see Kun ct al. 24) 



"From Kun et al. (2004). 

Covered only by our MIPS map. 
c Covered by our IRAC and MIPS maps. 

Has a mid-IR excess as determined here via Spitzer colors. 
°Lco ot al. J2^) identifies 2MASS 05073016-0610158 as an M0, and 2MASS 05073060-0610597 as an M2. 



Table 2. Summary of IRAC and MIPS observations 



Program 


AORKEY 


map center 


date obtained 


IRAC 


235 


13375232 


5:07:33.00 


-6:17:00.0 


2005-03-25 


266 


16940288 


5:05:52.00 


-6:27:52.0 


2006-03-24 


266 


16941056 


5:03:55.00 


-6:42:22.0 


2006-03-24 


266 


16941312 


5:05:17.00 


-7:10:38.0 


2006-03-23 


266 


16941568 


5:04:24.00 


-7:41:48.0 


2006-03-24 


266 


16941824 


5:03:10.00 


-7:04:57.0 


2006-03-24 


462 


24248832 


5:06:29.50 


-6:52:58.1 


2008-03-09 


462 


24249088 


5:07:18.39 


-6:07:11.1 


2007-10-16 


462 


24249344 


5:07:31.60 


-6:21:50.1 


2007-10-16 


MIPS 


235 


13374976 


5:07:33.00 


-6:17:00.0 


2005-03-10 


266 


16940544 


5:06:10.00 


-6:47.47.0 


2006-03-31 


266 


16940800 


5:03:42.00 


-7:29:19.0 


2006-03-31 



10 



by individual inspection. This resulted in a catalog of 84 point sources. The zero point used for 
conversion between flux densities and magnitudes for the 70 /xm data is 0.77 Jy, again from the 
MIPS Data Handbook. 



No point sources were detected in the 160 ^m data. 



2.4. Optical Observations and Data Reduction 

We obtained UVR C I C images by observing with the United States Naval Observatory (USNO) 
40-inch telescope (Flagstaff, AZ) over several epochs between November 2006 and January 2007. 
An area of 2.8 square degrees was covered in all 4 bands by mosaicking the ~23.9'x23.2' field of 
view. We acquired at least two images at two different exposure times for each pointing; see Table[3j 
On 15 March 2010, an additional epoch of relatively short (<60 sec) integrations was obtained in 
VR C I C just of the head of the nebula (where our YSO candidates are located; see below) as a final 
check on our measurements. Aperture photometry was performed on the objects in each frame, 
magnitudes between ~10-12 in the shallowest integrations and ~20-22 in each band in the deepest 
integrations, subject to variations in the quality of each night. Absolute photometric calibration 
was achieved through observation of Landolt standards during each night. 

The additional epoch in 2010 allowed us to check our calibration and assess intrinsic variability 
in these stars on ~3 year timescales. Several of the YSO candidates are variable, some highly 
variable (>0.2 mag; see Table 4). The three most variable sources are all embedded in reflection 
nebulae. 

In order to include the optical data in our spectral energy distributions, we needed to convert 
the measured magnitudes to flux densities in Janskys. For the U and V filters, we used the 
Johnson zero points from Allen's Astrophysical Quantities (Fourth Edition, 2001, Arthur N. Cox 
(ed.), Springer- Verlag) , specifically 1823 and 3781 Jy, respectively. For the R c and I c filters, we 
used the zero points from (J4J), which are 3040 and 2433 Jy, respectively. 



2.5. Bandmerging 

We matched IRAC sources from adjacent channels to empirically estimate our astrometric 
uncertainty. The histograms of coordinate differences show a lcr uncertainty of ~0".3 in both 
directions. At this Galactic latitude and longitude (207.4,-27.9), the density of sources is relatively 
low (~10 detections per square arcmin in the 3.6 /im images, 4 detections per square arcmin in R c 
and I c images). Via empirical tests across all bands, we determined that a 2" radius (1.6 IRAC 
pixels) for source matching works for our sources across the optical, 2MASS, IRAC, and MIPS 
bands. For a matching radius less than ~l", not all of the real matches are found. The number of 
matches increases steeply with radius until 1"; after that, the number of matches increases slowly 



Table 3. Summary of our UVR C I C observations from Nov. 2006-Jan. 2007 



# 


RA (2000) 


Dec (2000) 


U exp. (s) 


V cxp. 


(s) 


R c cxp 


00 


I c cxp. 


(s) 


U date 


V date 


Ft,. 


late 


Ic 


late 


1 


05:08:14.820 


-6:16:39.062 


1801 


, 180 


900, 


60 




20, 600 




60, 600, 


10 


18/11/06, 


17/11/06 


18/11/06, 


17/11/06 


17/11/06, 


18/11/06 


17/11/06, 
17/11/06 


18/11/06, 


2 


05:08:12.326 


-6:16:49.188 


900. 


180, 


900, 




ISO. 


5, 900, 


ISO. 


20, 30 


2, 


31/10/06, 


17/11/06, 


31/10/06, 


31/10/06, 


31/10/06, 


31/10/06, 


21/11/06, 


31/10/06, 








360 




60, 


120 


30 


20, 30, 


2 


5, 900, 
60, 600, 


ISO. 
10 


09/01/07 




20/11/06, 
21/11/06 


09/01/07, 


31/10/06, 
31/10/06, 


20/11/06, 
18/11/06 


18/11/06, 
31/10/06, 

20/11/06. 
21/11/06 


31/10/06, 
31/10/06, 
18/11/06, 


3 


05:08:11.814 


-6:04:28.112 


1801 


, 180 


900, 


00 




20, 600 




60, 600, 


10 


18/11/00. 


17/11/06 


18/11/06, 


17/11/06 


17/11/06, 


18/11/06 


17/11/06, 
17/11/06 


18/11/06, 


4 


05:06:54.960 


-6:13:32.816 


ISO. 
360 


600, 


900, 

121) 




60, 


20, 600 




60, 600, 


10 


17/11/06, 
09/01/07 


19/11/06, 


19/11/06, 
09/01/07 


17/11/06, 


17/11/06, 


19/11/06 


20/11/06, 

19/11/06 


19/11/06, 


5 


05:05:35.082 


-6:13:31.057 


ISO. 


600 


900, 


00 




20, 600 




60, 600, 


10 


17/11/06, 


19/11/06 


19/11/06, 


17/11/06 


17/11/06, 


19/11/06 


17/11/06, 
19/11/06 


19/11/06, 


6 


05:06:54.658 


-6:33:28.589 


ISO. 


(i(l() 


900, 


00 




20, 600 




60, 600, 


10 


17/11/06, 


20/11/06 


19/11/06, 


17/11/06 


17/11/06, 


19/11/06 


19/11/06, 
19/11/06 


19/11/06, 


7 


05:05:34.572 


-6:33:29.144 


ISO. 


(i(l() 


900, 


00 




20, 600, 


40 


20, 60, 
10 


iOo, 


17/11/06, 


21/11/06 


21/11/06, 


17/11/06 


17/11/06, 
09/01/07 


21/11/06, 


09/01/07, 
21/11/06, 


17/11/06, 
17/11/06 


8 


05:04:14.580 


-6:33:27.864 


180, 


000 


900, 


00 




20, 60, 


,00 


60, 600, 


10 


17/11/06, 


21/11/06 


21/11/06, 


17/11/06 


17/11/06, 
21/11/06 


21/11/06, 


21/11/06, 
17/11/06 


21/11/06, 


9 


05:07:25.289 


-6:48:31.295 


180, 


000 


1)00. 


00 




20, 600 




60, 600, 
2 


10. 


17/11/06, 


21/11/06 


21/11/06, 


17/11/06 


17/11/06, 


21/11/06 


21/11/06, 
17/11/06, 


21/11/06, 
21/11/06 


111 


05:06:05.416 


-6:48:28.150 


180, 
360 


600, 


1)00. 
120 




60, 


900, 

000 


20. 


15, 900, 
600, 10 


00. 


17/11/06, 
09/01/07 


21/11/06, 


21/11/06, 
09/01/07 


17/11/06, 


09/01/07, 
19/01/07 


17/11/06, 


09/01/07, 

17/11/06, 
17/11/06 


09/01/07, 
19/01/07, 


11 


05:04:47.396 


-6:48:48.174 


180, 


000 


1)00. 




60, 


900, 


20. 


900, 


00. 


17/11/06, 


19/01/07 


19/01/07, 


17/11/06, 


09/01/07, 


17/11/06, 


09/01/07, 


17/11/06, 












600, 


10 




600, 2 




600, 10 


2 






09/01/07, 


09/01/07 


19/01/07, 
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with radius until ~4" where the number of matches is dominated by false source associations. Using 
a matching radius of 2" appears to provide the best possible merging of sources across bands, even 
for 24 /j,m; the centroiding for the 24 /xm data is much better than 2", even though the data have 
a spatial resolution of ~6". 

Because the 3.6 /jm observations are the most sensitive of our Spitzer observations, we started 
by constructing a catalog based on these detections. To be included this catalog, sources have to 
be detected in 3.6 /jm and in at least one of the other bands. Out of the 39,817 sources in this 
catalog, 3365 (8.5%) have been detected in all 4 IRAC channels, 15% in the JHK S bands, 50% 
in the I c band (recall that the optical observations are deeper than the 2MASS observations), and 
just 12% in the U band. At 24 /iin, 1% (426 sources) of the 3.6 /mi-based catalog are detected. 
However, since the spatial coverage provided by our MIPS map is larger than that for IRAC, and 
since some very embedded young objects may not be detected in shorter wavelengths than 24 //m, 
we have also included in the catalog those sources with 24 (iva. detections even if they do not have 
3.6 /mi detections. Out of the 1082 sources detected in 24 /mi, 48% are detected in the I c band 
and 27% in the U band. However, the region mapped by MIPS is not exactly matched to the 
optical observations; considering just the region covered by both the optical and MIPS maps, these 
numbers become 69% and 39%, respectively. Similarly, 30% of the 24 /xm sources are detected in 
all 4 IRAC bands; if one is restricted to the area covered by both MIPS and IRAC maps, 70% are 
detected at all 4 IRAC bands plus MIPS-24. Finally, 38% of the 24 //m sources are detected in 
JHK q . 



3. Selection of YSO candidates with IR excess 

With our new multi- wavelength view of the molecular cloud, we can begin to look for the 
members of this complex. We focus on finding sources having an infrared excess characteristic of 
YSOs surrounded by a dusty disk. We do not expect to find a dense population of YSOs (see JTJ. 
Additionally, the molecular cloud is not very dense, so it will not obscure many of the background 
sources. Thus, we expect a high source contamination rate. Separating the YSO candidates 
from the extraneous sources requires an extensive weeding-out of the galactic and extragalactic 
contaminants. To begin this task, we identify objects with excesses at IRAC and MIPS wavelengths. 
Then, we use filtering mechanisms from the literature based on Spitzer colors to distinguish likely 
galaxies from likely members. We also use the additional information we have from the optical to 
further winnow the candidate list. 



3.1. IRAC selection 

We first select sources using a set of IRAC color-magnitude and color-color diagrams. Our 



selection is based on the (J20l ) method and is described in (|18l ): in summary, this method uses cuts in 
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several IRAC color-color and color-magnitude spaces to identify likely YSOs and likely contaminants 
such as active galactic nuclei (AGN). Using this technique, we identified 459 sources having colors 
consistent with galaxies dominated by polycyclic aromatic hydrocarbon (PAH) emission and 765 
more sources having colors consistent with AGN. Just 27 sources are not flagged as background 
contaminants and have colors compatible with YSOs with an IRAC excess in the [3.6] — [5.8]/[4.5] — 
[8] plane (where the bracket notation denotes magnitudes, e.g., [3.6] means magnitude at 3.6 /im). 
Figure |4] shows the IRAC color-color diagram for the sources that passed and failed these color 



tests. (We note for completeness that (|19l ) have updated their YSO selection method with new 



rejection criteria; we have verified that these new criteria do not affect our list of YSO candidates.) 

Since we expect to find a lot of contamination here, the list of 27 sources not flagged as con- 
taminants bears further scrutiny. These 27 sources consist of 9 sources brighter than [3.6]=12, with 
the remaining 18 all fainter than this limit. The SWIRE (Spitzer Wide-area Infrared Extragalactic 
Survey; Lonsdale et al. l3ll ) observations of the ELAIS Nl extragalactic fielco consists of integrations 



of similar depth, but this sample is expected to be essentially entirely galaxies and foreground stars, 
and as such provides a guide to the range of values expected from these contaminants. The SWIRE 
survey shows very few galaxies have [3.6] brighter than 12 magnitudes. The IC 2118 objects that 
we have identified explicitly as contaminants have 3.6 /im flux densities comparable to the SWIRE 
contaminants, with [3.6] ranging from 12 to 18 mags. But, out of the YSO candidates, 18 also have 
magnitudes fainter than [3.6] = 12. As the brightness decreases, the chances of these objects being 
background contaminants increases. 

Since we have a wealth of supporting data, we can use these data to additionally constrain our 
YSO selection. Figure [5] shows an optical color-magnitude diagram for these sources. The optical 
colors and magnitudes for these 18 faint sources all fall in the region occupied by main-sequence 
stars or background galaxies. In contrast, the 9 remaining YSO candidates with [3.6] < 12 all 
appear redder, near or above a 30 Myr Siess et al. (2000) isochrone in the V/V — I c color magnitude 
diagram. We strongly suspect that the faint ones are contaminants. The 9 brighter objects selected 
in [3.6] - [5.8]/[4.5] - [8] with [3.6] < 12 are included in our list of IRAC-selected YSO candidates. 
For reference, the 18 objects with [3.6] > 12 are represented with a blue filled circle in most of the 
remaining figures in this paper. The restriction we have imposed that [3.6] < 12 means that we are 
effectively senstitive to masses greater than 0.1 M & at a distance of 210 pc for an assumed age of 
~3-5 Myr. 

All of these 9 IRAC-selected YSO candidates are seen at MIPS-24. This is in contrast to the 
catalog as a whole, where only 10% of the sources seen in all four IRAC bands are also seen at 24 
/mi. Since legitimate members of IC 2118 with infrared excesses located at ~210 parsecs should be 
detected at 24 /im by our observations, the fact that we see these objects gives us further confidence 
that these objects are truly YSOs. Three of them are also seen at 70 /im; similarly, legitimate YSOs 
with large excesses should be detected by our shallow 70 /mi observations, so this also suggests that 



2 VizieR Online Data Catalog, 11/255 (|4 
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these objects are true YSOs. (Note, of course, that most of the sources seen in our observations at 
24 or 70 /jm are foreground stars or background galaxies, so detection at any MIPS band by itself 
is insufficient evidence for youth.) 

We extracted and examined the POSS, 2MASS, and our optical images for each of these 
candidates to verify that they did not appear extended in any of these bandsj. While all extended 
objects may not be galaxies (many YSOs in Taurus at 140 pc appear as extended objects due to 



circumstellar reflection nebulae, edge-on disks, jets, etc.; see Rebull et al. 1371 ). extended objects 
seen here are much more likely to be galaxies than YSOs at 210 pc. All of our YSO candidates 
passed this check and appear to be point sources in all available bands. 



3.2. MIPS selection 

This IRAC-based selection method does not allow us to detect young stars having an infrared 
excess only at wavelengths longer than 8 /zm. We looked for potential sources with a [3.6] — [24] color 
excess but no IRAC excess, which would be intermediate between photospheric objects and Class II 
stars; these objects are Class III sources with weak excesses (possible transitional disks). Figured] 
shows the [3.6]/[3.6] — [24] color-magnitude diagram. In this parameter space, photospheres have 
zero color, and galaxies are red and faint. This figure also contains data from the SWIRE ELAIS 
Nl field. We have truncated the SWIRE catalog to fit our shallower IC 2118 24 /xm detection 
limit, rejecting SWIRE sources with [24]>9.4, and plotted the catalog here as black contours. The 
SWIRE sample is essentially entirely galaxies and foreground stars, and as such provides a visual 
guide to the locations where such objects appear in the diagram. 

The most striking thing about this diagram for IC 2118 is that there are three distinct groups 
of objects: objects of zero color (likely foreground/background stars), objects that are faint and 
red (likely galaxies), and objects that are bright and red (likely YSOs). No sources appear between 
the photospheric/Class III objects and the Class II objects. The bright objects we selected above 
based on IRAC colors are also bright here, strongly suggesting the presence of an infrared excess 
we interpret as circumstellar dust. The contaminants identified above are also identified as likely 
contaminants here; these objects nearly all fall on top of the SWIRE contours. The objects identified 
as having YSO-like IRAC colors but being as faint as the contaminants are also faint and have 
contaminant-like colors in this parameter space. Within the faint contaminant portion of the 
diagram, the YSO candidates are among the brightest or reddest objects, which is a likely result of 
our selection process. This lends further support to our assertion that the faint objects are likely 
background galaxies. 



3 Some of our high school students, trained in observational techniques using the Sloan Digital Sky Survey Galaxy 
Zoo, participated in this classification. At least three different individuals examined images of each of more than 400 
candidates selected by a variety of Spitzer color selection mechanisms in order to identify likely galaxies. 
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[5.8] - [8] 



Fig. 4. — The [3.6] — [4.5] versus [5.8] — [8] color-color diagram. Sources nagged as PAH emission 
sources (background galaxies) or AGN by our selection method are plotted in green squares. Red 
star symbols are selected YSOs candidates with [3.6] < 12 and blue points are YSO candidates 
selected in the same fashion but with [3.6] > 12. The red triangle is an edge-on disk candidate 
(see M3.2|) . The black points are the remaining detections. The x symbol near zero color is the 
weak-lined T Tauri (wTTS) candidate (2MASS J05060574-0646151, see jJOJ). See text for more 
discussion of each of these categories. 
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Fig. 5. — Optical V versus V — I c color-magnitude diagram. The symbols are the same as in 
Figure [U Note that the blue dots representing the [3.6]>12 sources are in fact plotted, but they 
are in amongst the green squares and black dots well below the ZAMS; this we take to be further 
evidence that the objects represented by the blue dots are contaminants and not YSO candidates. 
The red triangle is the edge-on disk candidate, which appears here below the ZAMS, as appropriate 
for such an object. The solid blue lines are models from Siess et al. (2000) at 1, 10 and 30 Myr 
and scaled to 210 pc, where we have tuned the color-effective temperature relation such that the 
100 Myr isochrone matches that of the Pleiades single-star sequence ()45l : 
is also indicated. 
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A reddening vector 
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Despite the fact that Figure[6]does not reveal any new candidate YSOs near the IRAC-identified 
candidate YSOs, it does identify another source of particular interest. Figure [6] highlights a source 
(as a red triangle) that is red, with a [3.6] — [24] color of ~6.5 and a moderately faint [3.6]~12, which 
is among the brighter values for the likely contaminants. This source (named 050721.9-061152) has 
not strictly been selected as a YSO candidate or background contaminant by our IRAC-based 
method, but it is located in the border of our YSO selection area in the [3.6] — [5.8]/[4.5] — [8] 
diagram; it has a clear [3.6] — [5.8] excess but a low [4.5] — [8] excess. Somewhat surprisingly, 
it is clearly detected at 70 /mi. Since our 70 /mi survey is shallow enough that we would only 
detect IC 2118 YSOs if they have a large infrared excess, we take that detection, plus the marginal 
shorter-band color selections, as evidence that this object may very well be a YSO candidate, and 
we add it to our list. We discuss this object again below. 

We use [3.6J vs. [3.6] — [24] to look for YSOs, whereas many papers in the literature (e.g., 



Rebull et al. [37j and references therein) use the K s vs. K s — [24] color-magnitude space to find 



YSOs. In that diagram, the K s — [24] colors for M stars deviate significantly from zero ()15l ). so 
objects without disks but with a (possibly unknown) M spectral type will appear to have an infrared 
excess. The color departure from for M stars using [3.6] — [24] happens at later types than for 
K s — [24], which makes it a better choice to use for searches for small excesses in cases where the 
spectral type of the candidates is unknown. However, one can obviously only construct a [3.6] vs. 
[3.6] — [24] diagram for those regions where one has both IRAC and MIPS data. Since our maps 
do not match exactly, and since there are previously-identified T Tauris well off of the densest 
parts of the nebula, we also constructed a K s vs. K s — [24] color-magnitude diagram using 2MASS 
and MIPS data to see if any YSO candidates appeared off the edge of our IRAC maps. Two new 
sources appear as potential YSOs here, but by inspection of the images, both of these sources are 
extended, and one (2MASX J05015386-0745321) is a known galaxy. We do not include either of 
these objects on our YSO candidate list. 

One previously-identified T Tauri star (RXJ 0502.4-0744) is located in our MIPS map but not 
in the IRAC map, and is detected at 24 /mi; however this object does not show K s — [24] infrared 
excess. Its measured flux density is 1.3±0.12mJy at 24 /mi ([24]=9.36±0.1). 



4. Properties of the YSO candidates 

4.1. Magnitudes and Spectral Energy Distributions 

Coordinates and our measured magnitudes between U and 70 /mi for our 10 YSO candidates 
appear in Table |U Errors in the optical photometry for these 10 objects are generally dominated 
by the systematic errors, and as such are essentially all ~0.1 mag at U, 0.03-0.05 mag at V, R c , and 
I c , 0.02-0.03 mag at JHK S (where these values come from 2MASS), and 0.02 mag at all four IRAC 
bands. The 24 /zm errors are typically 0.04 mag; the 70 /mi errors vary enough between objects that 
they appear separately in the table. Spectral energy distributions (SEDs) for the 9 IRAC-selected 
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Fig. 6. — The [3.6]/ [3.6]- [24] color magnitude diagram. Symbols are the same as in Figured! The 
black curves are contours from the SWIRE catalog truncated to our IC 2118 detection limit (i.e., 
[24] < 9.4); see text. 
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YSO candidates appear in Figure [71 To guide the eye, we wished to add reddened stellar models to 
these plots, but spectral types are only known for four of the nine candidates. In order to provide 
a reference, for the remaining objects, we assumed an MO type. Thus, for each object in Figure [7J 



a reddened model is shown, selected from the Kurucz-Lejeune model grid (|29l ; l30l ) and normalized 
to H band. Note that this is not meant to be a robust fit to the object, but rather a representative 
stellar SED to guide the eye such that infrared excesses are immediately apparent. 

Our 10th YSO candidate is the MlPS-selected object, 050721.9-061152. This object's SED 
appears in Figure [8] and shows that the strong infrared excess starts after ~8 /im (which is why 
this object is not found via the IRAC-based selection method). This kind of SED is typical of 
a YSO seen edge-on. At wavelengths shorter than ~24 ^m, the thermal emission is absorbed by 
the disk itself and the source is mostly seen in scattered light; this object is indicated by the the 
triangles in FigureEl where it is below the pre-main-sequence locus defined by the other members. 
At longer wavelengths, the disk becomes optically thin and the thermal emission dominates the 
SED. In FigureJHl we have compared the SED of 050721.9-061152 with two well known edge-on disk 



SEDs: HH 30 (|8j) and HV Tau C (|13l ). The flux densities have been normalized to match that of 
050721.9-061152 in the J band. The figure illustrates that the global shape of the 050721.9-061152 
SED, with a minimum at 8 fj.ni, is compatible with an edge-on disk YSO. In our color-magnitude 
(Figures [5] [TOj and [6|) and color-color (Figures H] and [9]) diagrams, this source is plotted with a 
red triangle. Its location in the optical CMDs excludes its identification as simply a debris disk or 
transition disk. Since photospheric optical and NIR emissions are seen through the disk in scattered 
light, they are hard to predict and broad-band observations cannot tell us much about the nature 
of the central object. If it is an edge-on disk, then it joins a relatively short list of similar objects; 
such objects can tell us much about the physical properties of the occulting edge-on disk. Further 
study of this object is warranted. 

Out of our entire list of ten YSO candidates, six are new discoveries. Four of our candidates 



are already known as young members of IC 2118; they were discovered by ([22j; |24| ) as classical 
T Tauri stars (cTTS) using Ha objective prism plates and confirmed by spectroscopy to be young. 
They are 050711.6-061510, 050653.5-061712, 050730.6-061060 and 050730.2-061016. For one of 
these, 050711.6-061510, the SED appears 'disjoint' between the I c - and J-band observations; this is 
not uncommon in T Tauri stars, and is usually explained by the star having undergone significant 
variation in between the epochs of observations. In fact, this star is one of three that appear in 
our optical observations to be nebulous and highly variable (>0.2 mag) even within the timescale 
of our observations. 



In the spirit of the (I25l ) and (|5ll ) Class I/Flat/II/III system, we fit a line to all of the available 
data between 3.6 to 24 fim. Eight of the IRAC-bright YSO candidates are Class II T Tauri stars 
and one (050729.4-061637) is a flat disk. These classifications did not change if we instead fit just 
3.6 to 8 (jlvd.. We did not attempt a similar fit of the edge-on-disk candidate, since the object is 
seen through the disk, and the slope changes significantly depending on whether the MIPS points 
are included in the fit. 
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Fig. 7. — SEDs for 9 of our YSO candidates; the 10th is in Figure El The squares are our optical 
points, the circles are 2MASS, the triangles are IRAC, and the diamonds are MIPS. Error bars 
are indicated, and in most cases are smaller than the symbol. Spectral types are known for some 
objects from (]24l ). and those are indicated where possible. In cases where no type is known, an M0 
is assumed. Reddened models from the Kurucz-Lejeune model grid (|29i ; |3Cj) are shown for reference; 
see text. 
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Fig. 8. — SED of 050721.9-061152 (black filled symbols), compared to the well known edge-on disks 
HH 30 (solid line, circles; Burrows et al. |8|) and HV Tau C (dashed line, plus symbol; Duchene 
et al.ll3l). These SEDs have been normalized to 050721.9-061152 in the J band. For the black filled 
points, the symbols are as in the previous Figure - squares are our optical points, the circles are 
2MASS, the triangles are IRAC, the diamonds are MIPS, and the error bars are smaller than the 
points. 



Table 4. Measured magnitudes for Spitzer-selected YSO candidates in IC 2118. 



SST Name 


2MASS Name 


U a 


V 


Re 


Ic 


J a 


H <T 


K s CT 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


[70] 


050650.1-061908 


05065009-0619079 


... ± ... 


21.59 


19.35 


17.43 


13.30±0.027 


12.16±0.026 


11.53±0.025 


10.75 


10.38 


10.04 


9.52 


6.79 


... ± ... 


050653. 5-061712 b 


05065349-0617123 


20.99±0.102 


17.00 


15.20 


13.71 


11.16±0.026 


9.91±0.027 


9.13±0.024 


8.20 


7.76 


7.31 


6.51 


3.27 


0.53±0.044 


050711. 6-061510 b '' 


: 05071157-0615098 


18.91±0.033 


16.29 


15.01 


13.85 


13.02±0.026 


11.28±0.027 


10.02±0.029 


8.74 


8.11 


7.63 


7.07 


4.78 


... ± ... 


050720.6-061612 


05072058-0616119 


21.93±0.225 


19.50 


17.58 


15.47 


12.91±0.020 


12.07±0.021 


11.63±0.021 


10.82 


10.30 


9.86 


9.13 


6.42 


... ± ... 


050721.9-061152 


05072188-0611524 


... ± ... 


19.36 


18.04 


16.88 


15.04±0.043 


13.99±0.045 


13.28±0.030 


12.27 


11.67 


11.36 


11.15 


5.45 


0.80±0.058 


050728.7-061326 


05072865-0613258 


20.29±0.058 


18.32 


16.68 


14.81 


12.50±0.028 


11.52±0.030 


11.04±0.024 


10.57 


10.17 


9.87 


9.25 


6.19 


... ± ... 


050729.4-061637 


05072937-0616368 


21.43±0.152 


19.03 


17.55 


16.10 


13.71±0.023 


12.41±0.024 


11.54±0.019 


10.27 


9.44 


8.80 


7.63 


4.33 


. . . ± ... 


050730. 2-061016 b '' 


= 05073016-0610158 


17.10±0.030 


14.99 


13.72 


12.54 


10.82±0.023 


9.57±0.027 


8.61±0.019 


7.30 


6.77 


6.35 


5.55 


2.43 


-0.39±0.030 


050730. 6-061060 b '' 


= 05073060-0610597 


17.58±0.030 


14.74 


13.47 


12.21 


10.12±0.025 


9.04±0.029 


8.25±0.015 


7.40 


7.02 


6.70 


6.13 


3.41 


1.20±0.103 


050733.7-061348 


05073372-0613474 


20.74±0.080 


19.19 


17.41 


15.39 


13.01±0.020 


12.39±0.026 


12.02±0.019 


11.40 


10.92 


10.55 


9.71 


7.08 


... ± ... 



a Thc errors in the VRI C1 IRAC, and MIPS-24 data arc dominated by systematic errors; they arc 0.03-0.05, 0.05, and 0.04 mag, respectively (sec text). 

b Prcviously-identified YSO; see Table[T] 

"Nebulous and highly variable (>0.2 mag) in the optical in our observations. 
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Fig. 9. — J — H versus H — K s diagram for the sample, with the same notation as earlier figures. 
The T Tauri locus is also indicated (dash-dot line) , as is a reddening vector (the arrow) . The main 
sequence is drawn by a solid line and the giant branch by small black triangles. The x symbol is 
the wTTS candidate (2MASS J05060574-0646151, see flO) . This diagram shows that all but one 
of the YSOs have an infrared excess starting at H-band and that they have moderate reddening. 
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Since spectral types are not known for six of the ten YSO candidates, it is not easy to determine 
masses, ages, or degree of reddening. However, Figure [9] shows the J — H versus H — K s diagram 
for the sample, with the same notation as earlier figures. By comparison to the ZAMS, it can be 
seen that all have moderate reddening. The source that does not show significant IR excess in the 
J — H versus H — K s diagram is located on the 1 Myr model isochrone. Likewise in the H versus 
H — K s diagram, Figure flOl the T Tauri stars appear well separated from the main sequence stars 
and the identified background contaminants. They all appear redder than the 1 Myr isochrone 
model mostly because they have infrared excess starting at -ff-band. If one makes the simplifying 
and probably incorrect assumption that the offset is entirely due to A v , and trace the objects back 
to the 1 or 3 Myr isochrone, the brightest objects are ~2M , and the faintest are ~0.1M Q . 



4.2. Previously-identified YSOs not recovered by our selection 



We note that (|24l ) reported other classical T Tauri star (cTTS) candidates in this region (see 
Table [1]) , but they are either on the southern end of the molecular cloud or significantly off the 
nebular region, and are not covered by our survey. 



d2J) reported a weak line T Tauri star (wTTS) candidate (2MASS J05060574-0646151) that 
may or may not be a true member of IC 2118. This object is within our surveyed region. This 
object does not show any infrared excess in the IRAC bands and is not detected at any of the MIPS 
bands. Therefore, our selection methods based on infrared excesses will not recover it. However, 
its position in our R c /R c — I c and V/V — I c CMDs is compatible with a pre-main sequence star 
according to models and the position of our 9 (IRAC-selected) YSOs in the same diagram. This 
object is represented with a x symbol in optical and NIR color-color and color-magnitude diagrams 
(Figures EJ [9j [TO]) and in Figure HU 



4.3. Location of the YSO candidates 

All 10 of our YSO candidates, including the edge-on-disk candidate, are located in the head of 



the nebula, i.e., G206. 4-26.0, the most massive molecular cloud (|22j ) of IC 2118. Figure [TT] shows 
the spatial distribution of the YSOs, using the same notation as in earlier figures. The distributions 
seen here lend further support to our assertions above that the green points are the contaminants 
(which are evenly distributed), the blue points (faint IRAC-selected objects) are likely galaxies 
(generally off-cloud and evenly distributed), and the IRAC- and MlPS-selected bright objects are 
likely YSOs. We expected to find YSO candidates further South in this cloud; however, evidently 
the conditions do not support substantial star formation in regions other than the "head." The 
surface brightness of the ISM at 8, 24, 70, and 160 //m is fainter in the southern portion of the 



nebula compared to the more northern portion. The CO maps from Kun et al. ()22j ) find that the 



head of the nebula, G206. 4-26.0, at 85 M Q , is at least 3 times more massive than any of the other 
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Fig. 10.— H versus H — K s diagram for the sample, with the same notation as earlier figures. 
The solid blue lines are models from (|42j ) at 1, 10 and 30 Myr and scaled to 210 pc. The x symbol 
is the wTTS candidate (2MASS J05060574-0646151, see $Oj) . 
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clouds. G206.8-26.5, with 28 M , houses only 2MASS J05060574-0646151, the wTTS candidate 
mentioned in the prior section, which may or may not be a true member of the cloud and does not 
have an infrared excess. At 24 ^m, the peak surface brightness of G206. 4-26.0 is ~34 MJy sr _1 , and 
the peak surface brightness of G206.8-26.5 is not that different at ~33 MJy sr _1 . At 160 //m, the 
peak surface brightnesses are ~205 and ~125 MJy sr _1 , respectively. The remaining clouds within 
our map, G207.3-26.5, G207.2-27.1, and G208.1-27.5, are 2, 2, and 6 M Q , respectively. Kun et al. 



(J22J) finds them all to be fainter, cooler, and physically smaller on the sky than either of the other 
two clouds in the head of the nebula. At 24 /mi, their peak surface brightnesses are ~31, ~31, and 
~27 MJy sr _1 , respectively; at 160 /xm, the peaks are ~77, ~30, ~80 MJy sr _1 , respectively. All 
of this suggests that the local conditions are different enough to significantly affect star formation 
in the two locations. 



4.4. Ultraviolet excesses 



Ultraviolet (UV) excesses can also be an indication of youth (see, e.g., Rebull et al. [38| and 
references therein). We have [/-band measurements for 9 of the known plus new candidate IC 2118 
members, including the WTTS; the two objects without U measurements are 050721.9-061152 and 
050650.1-061908 (see Table H]). Because these objects are subject to moderate reddening, which 
will significantly affect the [/-band observations, a robust determination of the UV excess requires a 
spectral type and secure determination of A v . However, the photospheres plotted to guide the eye in 
Figure [7] suggest that there are UV excesses in at least 6 of our 9 IRAC-selected YSO candidates. 
Two of the largest apparent UV excesses are in two of the stars with known types, 050730.2- 
061016 and 050730.6-061060 (omitting 050711.6-061510 because of the evident stellar variation), 
suggesting that, when better spectral types are obtained for the remaining objects, UV excesses 
may be apparent in most of them. 



4.5. On the Distance to IC 2118 

At least two (perhaps) intertwined mysteries enshroud the Witch Head - its distance and the 
external source responsible for sculpting its surface, illuminating it, and possibly triggering star- 
formation within it. There is at least a factor of two uncertainty in the distance to IC 2118 - 210 
pc vs. 440 pc. Our process for identifying YSOs towards IC 2118 does not depend strongly on 
the distance; the same stars would have been selected for any distance up to the point where the 
faintest objects start to blend with the background galaxies, and all of our YSOs are considerably 
brighter than this limit. However, physical interpretation of the IC 2118 YSOs depends on accurate 
determination of the distance. For that reason, we return to the distance question now. 

Two sources (at different distances) are listed as potential illumination and/or wind-sculpting 
sources for the nebula - the Trapezium and Rigel. Figure 3 includes two arcs from circles centered 
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Fig. 11. — Distribution of our YSO candidates and background galaxies or AGN. Symbols are 
identical to those in Figure gj The x symbol is the wTTS candidate (2MASS J05060574-0646151, 
see §4.2p . The background gray scale image is the 160 /im MIPS observation of IC 2118. 
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Fig. 12. — A 3-color composite of IRAS data at 12, 25, and 100 //m, centered on IC 2118, with 
a circle at the location of 9 l Ori C and Rigel. The region shown is ~19° on a side. This figure 
illustrates that IC 2118 is one of several nebulae located many degrees from the ONC whose 
morphology appears "wind-blown" by a source roughly towards the ONC and belt stars. 



29 




V-l 



Fig. 13. — Optical My versus V — I c color-magnitude diagram. The (|42l ) isochrones are included 
(1, 10, and 30 Myr), but shifted to absolute My. The black stars/triangle are our YSO candidates, 
same notation as earlier figures, assuming a distance of 440 pc, and the red stars/triangle are our 
YSO candidates, assuming a distance of 210 pc. The grey x symbols are Taurus YSOs (from Rebull 
et al. 2010 and Giidel et al. 2007 and references therein), taken to be at 140 pc; Torres et al. 1471.148. 
The Taurus distribution is broad and there are many fewer IC 2118 stars, but this distribution 
weakly suggests that IC 2118 is closer rather than farther (see text). 
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on the Trapezium (7.35° from IC 2118) and Rigel (4.48° from IC 2118). We discuss these more 
below. 



4-5.1. Arguments in Favor of d ~ 400 pc 
The most secure and accurate distance for any part of the Orion complex is the distance derived 



using the VLBA for several members of the Orion Nebula Cluster (ONC) by (|34l ). which is d = 414 
± 7 pc. As the most massive members of the ONC, we assume that the Trapezium cluster stars, 
and in particular 8 l Ori C, are at this distance. Presumably, this also means that the Orion Nebula 
is at this distance. At slightly lower accuracy, the BN/KL complex just north of the Trapezium 
has been determined to be at 437 ± 19 pc via VLBI measurement of water masers. This puts the 
most active current star formation and the densest molecular gas in Orion at > 400 pc. 

As noted by (|2j) and expanded upon considerably by Ogura & Sugitani (35j,hereafter OS98), 



there is a large population of nebular structures surrounding the ONC, many of which have cometary 
shapes or other morphology suggestive of shaping by winds or UV photoionization. If arrows are 
attached to each structure pointing towards the apparent direction of the wind, these arrows form 
a centro-symmetric pattern pointing to an origin roughly between the Trapezium and the Orion 
belt stars (see Figure 2 of Bally and Figure 2 of OS98). Both Bally and OS98 include IC 2118 as 
a nebula that exemplifies this pattern. Figure [T2"| based on IRAS data, shows a wide-angle view of 
the gas and dust west of Orion (centered on IC 2118), illustrating that IC 2118 is simply one of a 
number of nebulae located many degrees from the ONC whose morphology appears "wind-blown" 
by a source roughly towards the ONC and belt stars. Detailed examination of the leading edge 
of the MIPS and IRAC images of IC 2118 (Figures [2] and [3]) supports this assertion. The arcs in 
Figure [3] show that while the morphology seems to better follow the arc from Rigel, the direction 
of the "windblown" structures in the nebula suggests that the Trapezium is instead the exciting 
source, or perhaps a structure centered near it, such as the expanding Orion-Eridanus Bubble (see, 
e.g., Kun et al. 2001). These arguments therefore support a distance to IC 2118 of order 400 pc. 



4-5.2. Arguments in Favor of d ~ 210 pc 

Rigel (V ~ 0.3, spectral type B8I) is located only about 2.5° from IC 2118, compared to the 
Trapezium and belt stars at more than 7° distant. In the optical, IC 2118 is normally described as 
a reflection nebula due to its blue colors and lack of a bright rim. The most obvious photon source 
- and the one most frequently mentioned in the literature - for illuminating IC 2118 is Rigel. Since, 
according to Hipparcos, Rigel is at d ~ 264 ± 25 pc, it follows that IC 2118 should be at a distance 
of order 250 pc. 



As noted by (|22l ). HD 32925 is located along the line of sight to IC 2118 and appears to have 



significant reddening (Ay = 1.36). This suggests it lies within or behind IC 2118. If one adopts 
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the Walraven photometric distance for HD 32925 (|7|) of d ~ 260 pc, this places IC 2118 nearer than 
this. However, the distance to HD 32925 is not well-established. The best spectral type available 
for HD 32925 is AOIII (Houk, Michigan Spectral Catalog); given V = 9.2, My = 0.0 Q and the 
Ay, this yields d ~ 370 pc. The original Hipparcos catalog parallax for HD 32925 was it = 1.69 ± 
1.35 mas, hence supporting any distance more than about 250 pc. The new Hipparcos reduction 
yields ir = 3.66 ± 1.2 mas, or d = 273 (+133, —66) pc. No firm conclusion can be drawn from 
these disparate estimates. 

Figure [131 presents the My vs. V — I CMD, comparing IC 2118 stars to YSOs in Taurus (with 
data from Rebull et al. 2010, Giidel et al. 2007, and references therein). Based on morphological 
grounds (e.g., the degree to which the YSOs are still embedded in their natal gas), we expect that 
the IC 2118 stars might be slightly younger than the often more physically dispersed Taurus stars. 
On the other hand, based on the ratio of Class I to Class II sources, the IC 2118 objects might be 
slightly older than Taurus. In Figure [T3l the distribution of IC 2118 objects is broad, but assuming 
a distance of 440 pc, then IC 2118 appears to be slightly younger than Taurus. Assming a distance 
of 220 pc, the bulk of the IC 2118 distribution is comparable to the bulk of the Taurus distribution 
at ~3 Myr. Figure [T3l thus weakly suggests that IC 2118 is closer rather than farther. A similar 
calculation with ONC members yields similar results. 



4-5.3. No clear answer! 

It is tempting to try to use the infrared colors to set some constraint on the distance, since 
after all the infrared is reprocessing the UV starlight, and the distance between a given source and 
the cloud clearly affects the radiation field bathing the cloud. The integrated (over the entire cloud) 
infrared spectral energy distribution is shown in Figure HU where for each band, a corresponding 
background outside the cloud has been subtracted and the the images have been smoothed to the 
lowest angular resolution before the measurements. The figure includes a fit performed using the 
DUSTEM model using the interstellar standard radiation field (ISRF; Mathis et al. 1983) In a 
nutshell, the DUSTEM model requires at least three different dust populations (PAHs, Very Small 
and Big Grains) to interpret the infrared dust emissivity (for details see Compiegne et al. 2008, 
2010; Desert et al. 1990); plus the incident radiation field that is the sum of the mean ISRF 
and a stellar component diluted by the distance squared. The fact that the observations are well- 
matched by the standard radiation field immediately indicates that neither 9 Ori nor Rigel are 
dominant sources in illuminating the nebula, and therefore, no simple contraint arises from the 
mid-IR measurements. 

Unfortunately, it is our opinion that the evidence currently available pertaining to the distance 
to IC 2118 is not definitive. GAIA could provide a definitive answer by deriving accurate parallaxes 
for HD 32925 and/or for the brighter YSOs of IC 2118. A spectral type for the diffuse optical 
emission from IC 2118 could establish whether the dominant photon source is Rigel or the hotter, 
earlier type stars of the ONC or belt region. 
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Fig. 14. — Best model fit to the extended dust emission infrared SED of the IC 2118 cloud itself; 
see text. The observed SED is plotted with circles and error bars, the modeled SED with boxes. 
Both are normalized to 160 /im. The solid line is the total dust emission spectrum; the other 
five lines are the different dust components (PAHs [dotted and medium dashed lines with emission 
features, from the cations and neutrals], stochastically heated Very Small Grains [long-dashed line, 
warmest blackbody curve], Big Grains at thermal equilibrium [long-dashed cool blackbody curve]; 
the smooth curve that peaks around 2-3 /mi [short dashed line] is a component added in to explain 
the near-IR continuum [see Flagey et al. 2006 and references therein]). The observations are well- 
matched by the standard radiation field, suggesting that neither 6 Ori nor Rigel are dominant 
sources in illuminating IC 2118. 



33 



4.6. J05080318-0617141: a high proper motion star 

While constructing 3-color images using Spitzer and POSS images, we noticed one object 
projected onto the nebulosity that appeared to have moved significantly between the epoch of the 
POSS image and the current epoch; see Table This object, J05080318-0617141, moves at a 
steady rate of ~0.2 arcseconds yr _1 between the 1955 and 2000 epochs; at a distance of 210 pc, 
this would be ~180 km s _1 . The object does not appear to move much between the 2MASS and 
IRAC epochs (2000-2005); five years may not have been enough time to see significant motion. 
Our optical observations (from 2007) do not have precise enough astrometry for this calculation. 
However, this object is in the USNO-B catalog of objects with measurable proper motion; it is 
object 0837-0052579, having /j, a =134 mas yr _1 and fj,g = —86.0 mas yr" 1 , quite comparable with 
the net ~0.2 arcseconds yr _1 noted above. 

We were granted two nights at the Palomar 200" in Jan 2007 to obtain follow-up spectroscopy 
of our YSO candidates. Unfortunately, due to high velocity winds, high clouds, and poor seeing 
on both of our nights, little of our data were usable. However, the spectrum we obtained of this 
relatively bright, high proper motion object reveals it to be an M star. 

We do not find an IRAC excess for this object, and it is not detected at 24 //m. We have optical 
and near-infrared measurements: [7=18.77, U=15.98, R C =U.62, I c =13.19, J=11.68, #=11.08, and 
-RT S =10.85 (where the near-infrared points come from 2MASS, and the error on all the optical points 



is ~0.03 mags). Using the V/V — I c diagram and the models from (1421 ), it is likely to be an M2 
to M5 star between 20 and 70 pc away. Assuming those distances, its projected space motion is 
between ~20 and 60 km s^ 1 . 



5. Discussion and Conclusions 

By using Spitzer to search for objects with infrared excess in the IC 2118 region, we have 
doubled the inventory of YSOs and YSO candidates, in the region for which we have both IRAC 
and MIPS data, from 5 (4 CTTS and 1 WTTS) to 11 (6 new candidates selected via Spitzer colors, 
plus the 4 CTTS and 1 WTTS from before). Including the YSOs found by Kun and collaborators 
outside our surveyed region, the total number of YSOs found in this region is between 15 and 17, 

Table 5. Observations of J05080318-0617141 



Data set 


position (J2000) 


date of observation 


(UT) 


notes 


POSS-1 


05 


08 


02.77-06:17:10.0 


1955-12-15 07:50 






POSS-2 


05 


08 


03.00 -06:17:11.9 


1985-01-20 11:26 






2MASS 


05 


OS 


03.19 -06:17:14.1 


2000-10-12 09:04 




JHK S images have a nearby "persistence artifact" 


IRAC 


05 


08 


03.18 -06:17:14.1 


2005-03-25 00:02 




large PSF 


our optical 


05 


08 


03.2 -06:17:14 


2007-01-09 02:50 




large errors on astrometry 
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depending on the membership status of two objects from (|24l ) whose membership is questionable 



Most of the confirmed YSOs, and all of the new objects, are found in the "head" of the nebula 



the cloud G 206.4—26.0 from (J22i ). This cloud is estimated to be 85 M Q total molecular mass (J22l ). 
If we assume that all 10 YSOs located in the head of the Nebula are legitimate members of the 
cluster and on average 0.5 Mq, then ~5% of the cloud's mass is in stars, which is a slightly higher 
star formation efficiency than could have been calculated before. 

The SEDs of all 10 of these objects are similar (Class II); 050729.4-061637 is the exception, 
with a flat SED that is much different than the others. The edge-on disk candidate 050721.9-061152 
has only a reasonably small excess until 24 fxm. All 10 of these objects are within ~13'; assuming 
a distance of 210-400 pc, these objects are all within ~0.8-1.6 pc of each other. 



In Perseus and other regions (Rebull et al. [39| and references therein), several small groups 



of likely young stars are found within a parsec of each other. Notably, they have a large range 
of SED types, in one case (Per 6) ranging from apparent photospheres to Class candidates and 
an apparently starless millimeter continuum core. The fact that the range of SED types is much 
smaller in the case of IC 2118 may simply be a result of the much lower extinction towards IC 2118, 
or that the stars in IC 2118 are old enough that any differences in the timescales of very early 
star formation (e.g., as a result of different initial disk masses or rotation rates) are wiped out, or 
it could be telling us something about the onset of star formation in these two different regions. 
There is no clear evidence that the stars in Perseus are a result of triggered star formation, but the 
literature suggests that IC 2118 was triggered. A wide range of SED types can also be found in 
the cometary cloud L1251 (see, e.g., Lee et al. 2006), suggesting that perhaps triggering may not 
result in a small range of SED types. 

If this region is indeed triggered star formation, then we might be able to find trends of age 



or mass with location. For example, ([24]) asserts that stars closer to center of the core should be 
younger. Since there are so few objects, since the distance is uncertain, and since we do not have 
spectral types for most of these objects, we can not do this rigorously; we can look at size of infrared 
excess, which is often taken as a proxy for age, e.g., in the Class 0/I/flat/II/III scheme. There is 
no clear correlation of various measures of IR excess and location in the cloud. This must be, at 
least in part, because so many of the SEDs are similar. 

Additional follow-up spectroscopic data are needed in the IC 2118 region to confirm or refute 
the YSO status of the six new objects with infrared excesses, some of which appear to also have 
ultraviolet excesses, and all of which have optical magnitudes compatible with being YSOs in 
IC 2118. The new edge-on disk candidate in particular warrants further study, since such objects 
are relatively rare. 

The distance to IC 2118 is still uncertain, with evidence to be found for both ~210 and ~440 
pc. GAIA will be needed to resolve this mystery. 
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